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Abstract 

It has been suggested recently by Gan, Coleman, and Andrei that 
studying the underscreened Kondo problem may help to understand 
the nature of magnetism in heavy fermion systems. Motivated by Do- 
niach's work on the 5 = 1/2 Kondo necklace, we introduce the un- 
derscreened Kondo necklace models with S > 1/2. The underscreened 
Kondo necklace is the simplest lattice model on which the competi- 
tion between Kondo spin compensation, and magnetic ordering due to 
an RKKY-type interaction can be examined. We used the mean-field 
approximation to determine the phase diagram, and found that the low- 
temperature phase is always an x — y antiferromagnet. This contention 
is further supported by the derivation of the exact form of the effective 
hamiltonian in the limit of very large Kondo coupling: it is found to be 
an antiferromagnetic x — y model for the residual S — 1/2 spins. In gen- 
eral, the degree of moment compensation depends on both the Kondo 
coupling, and on S. 
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1 Introduction 



The description of the magnetically ordered states of heavy fermionic systems 
is a problem of great current interest Q. It is usual to find that the ordered 
moment is rather smaller, and in quite a few cases much smaller, than what 
would correspond to isolated 4/ (or 5/) shells. On the other hand, at elevated 
temperatures, the paramagnetic susceptibility is the one expected for isolated 
/-shells in a crystal field. The phenomenon of reduced moment magnetism 
at low temperatures is supposed to result from the competition of the Kondo 
effect, and RKKY interaction 0. The Kondo effect strives to set up a state in 
which the spins of /-electrons are compensated by those of conduction elec- 
trons; the RKKY interaction tends to order the localized moments, preventing 
the formation of an overall singlet state. The outcome of this competition can 
be apparently anything between the extremes of a conventional RKKY magnet 
with well-defined localized moments, and a non-magnetic heavy Fermi liquid. 
We can envisage this by saying that the Kondo effect progresses up to a point, 
reducing the moments to a fraction of the ionic value, and then the residual 
moments get ordered. - - The description of this situation is greatly compli- 
cated by the fact that the Kondo compensation clouds corresponding to nearby 
ions strongly overlap [[| ffl : there is a continuing controversy as to whether the 
collective Kondo effect of the lattice has an energy scale essentially different 
from the single ion Kondo energy. 

The treatment of the competition between magnetic and non-magnetic 
states is apparently quite difficult in the case of models with the the feature of 
perfect screening: 25* = k, where k is the number of screening channels, and 
S is the localized spin. As demonstrated in detail M for the case S = 1/2, 
this perfect balance allows the construction of a strictly non-magnetic heavy 
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fermi liquid; it remains an open question whether such a non-magnetic state 
is necessarily unstable against magnetism [pfl and/or superconductivity 0. 

It has been suggested recently by Gan, Coleman, and Andrei [[/J that one 
should be able to gain an insight into the nature of heavy fermion magnetism 
by studying the underscreened Kondo problem with 2S > k. The single- 
ion problem is solved exactly ||, demonstrating that in the ground state, a 
partially compensated spin S — k/2 remains. The underscreened Kondo lattice 
(in standard notations) 

H KL = £ e k c+ Q c fcQ + ^ E E e^-^cL**^ ■ Si (1) 

ka a,/3 kk'j 

with an orbitally non-degenerate conduction band (k = 1), and S > 1/2, seems 
destined to become magnetic (near half-filling, presumably antiferromagnetic), 
but nothing specific is known yet. 

Recalling that in the perfectly screened case k = 1, S = 1/2, the simplest 
solution of the Kondo necklace model introduced by Doniach || has been a 
useful guide to studying the more complicated problem of the Kondo lattice, 
we suggest that the introduction of the underscreened Kondo necklace models 

H KN = jJ2S y r i + W EOfrf + if if) (2) 
j (id) 

should be helpful as a prelude to a study of the underscreened Kondo lattice. 

The great simplification achieved by replacing the Kondo lattice model (1) 
by the Kondo necklace model (2) is due to the fact that H^l has in it spins 
and fermions, while Hkn is formulated solely in terms of spin operators. The 
spins S are the same localized spins as in (1). We consider an arbitrary spin 
S; the underscreened models belong to S > 1/2, but to facilitate comparison 
with the case of exact Kondo screening, we include also results for S — 1/2 
which belongs to the original necklace model. 
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The pseudospins r, with |r| = 1/2, are meant to represent the spin degrees 
of freedom of the conduction electrons. Charge degrees of freedom are omitted 
on the ground that they are likely to belong to higher energies, and are thus 
considered to be of little relevance for the low-energy world of Kondo physics. 
Note that one r-spin for each site should really correspond to a half-filled 
conduction band. At this special filling, the ground state of Hkl is expected to 
be insulating: either because of Luttinger's theorem, in case it is non-magnetic; 
or because of magnetic cell doubling, in case it is an antiferromagnet. In either 
case, there is a charge gap, providing a further justification for omitting the 
charge degrees of freedom in H KN . 

In one dimension (Id), the x — y model of the r-spins can be mapped by 
the Jordan-Wigner transformation to a tight-binding model of spinless non- 
interacting fermions. In this case, it becomes quite apparent that the x — y 
term in (2) stands for a sea of fermions. However, even in higher dimensions, 
it is clear that Hkn incorporates the competing tendencies of magnetic order- 
ing, and local Kondo singlet formation. This being just the question we are 
interested in, we do not hesitate to consider the model (2) for lattices of arbi- 
trary dimensionality and structure (even though the term "necklace" becomes 
a misnomer). — In any case, the simple mean-field treatment which we carry 
out here, in close parallel with Doniach's first study || of the original necklace 
model, is really justified only for 3d, and gives at best a qualitative guide as 
to what to expect in lower dimensions. 

The first term in (2), with the antiferromagnetic Kondo coupling J > 0, 
favours the formation of local low-spin states with \S + t\ = S — 1/2. The 
second term, with W measuring, loosely speaking, the kinetic energy of the spin 
degrees of freedom of the conduction electron sea, couples nearest-neighbour 
sites, and thereby mixes in also components from the local \S + t\ = S + 1/2 
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high-spin states. 

In the perfectly screened case S = 1/2, Doniach || found a ground state 
phase transition between an antiferromagnetic state at J < J cr = W, and a 
fully Kondo-compensated overall singlet state at J > W. The antiferromagnet 
has gapless spin excitations, while the Kondo state shows a spin gap. — It is, 
of course, always questionable whether mean-field results can be trusted, and 
it has been a matter of long debate what the true behaviour should be like. The 
problem is most subtle in Id, for the true "necklace". Amazingly, the mean- 
field approximation turns out to be a good guide even in this case. Naturally, 
true long-range order has to be replaced with quasi-long-range order, but it 
seems to be true that there is a ground state phase transition between a low-J 
state which is almost magnetic, with the spin excitation spectrum gapless, and 
a high- J Kondo state where there exists a spin gap. In finding this, very recent 



numerical studies JT0[] corroborate earlier renormalization group results |p"Tf , 
refuting some Monte Carlo work [|TIJ which claimed that the spectrum is gapful 
for all J > 0. Note, however, that for the Id Kondo lattice, it is reasonably 
expected, and also convincingly demonstrated |13] , that the spectrum is always 
gapful. This indicates that the breaking of the full spin-rotational invariance, 
which is caused by introducing the x — y form of the intersite coupling in Hkn, 
puts the Kondo necklace into a different universality class than the Kondo 
lattice. — In the present context of the underscreened Kondo necklace model, 
we will understand this more clearly after deriving the effective Hamiltonian 
governing the large-J behaviour which turns out to be an antiferromagnetic 
x — y model of the partially Kondo-screened residual moments. 
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2 X— Y Type Ant iferro magnetism in the Neck- 
lace Model 



Having thus found sufficient justification for doing simple mean-field theory to 
get a first impression of the behaviour of the model (2), we proceed to decouple 
the x — y term, by assuming a <i-dimensional cubic lattice, and postulating a 
two-sublattice (A and B) antiferromagnetic order of the r-spins 

/ t > if i E A 
{Ti) = { -t if i E B (3) 

while (r y ) = everywhere. 

The single-site mean-field hamiltonian is 

h MF = J St - 2ut x (4) 

where uo = Wtd, and t = (r x ) is the thermal expectation value of r x , which 
has to be determined self-consistently from 

E |£f e _ e , /T W 

Here €j, and are the eigenvalues, and eigenstates of h MF . 

Conveniently chosing the quantization axis as the rr-axis, the 2(2S + 1)- 
dimensional eigenvalue problem immediately separates into 2S two-dimensional, 
and 2 one-dimensional problems. In the subspace S x + r x = S — m + 1/2, 
where m — 1, 2, ..^S*, the matrix of h MF is found to be 

' j(S-m)-uj ±y/m(2S -m+1) N 

v 1^(25-771+1) -^(5 -771+ 1) +U) J 

There is no point to writing down the solution of the diagonalization, we pro- 
ceed straight to the results obtained from solving the self-consistency equation. 



(6) 
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The Neel temperature Tjv is obtained from solving the linearized self- 
consistency equation 



Wd r 1 
3Z {1 + 2S)^T~ 



N 



5(25 



J(S+1) 

l)e~^~ + (5+ l)(25 + 3)e~ 



JS 
2T N 



+ 



1 165(5+1) 



J 1 + 25 



J(S+1) 



where 



J(S+1) 

5e ^ + (5 + l) e " 



2T N 



(7) 



(8) 



Some representative curves giving the antiferromagnetic-paramagnetic phase 
boundary in the J/W — T^/W plane are shown in Fig. 1. At J=0, all curves 
start from the common value T/v = Wd/2, the mean-field solution for the x—y 
model of the r-spins. This single point is, of course, somewhat pathological, 
since the decoupled 5-spins remain completely disordered at all temperatures, 
including T = 0. For all J > 0, however, the antiferromagnetism of the r-spins 
induces an (oppositely polarized) antiferromagnetic order of the S-spins. 

As the Kondo coupling reaches the range J ~ Wd/S, the states derived 
from | S + t\ = 5 + 1/2 levels of the single-ion problem are beginning to make 
a negligible contribution to thermal expectation values, and T N quickly drops 
toward its asymptotic value 



Wd 25-1 



lim T N = (9) 

J/Wr-oo 6 25+1 K J 

For 5 = 1/2, there is no antiferromagnetism at large J's, in agreement with 
the solution found by Doniach |J. However, for all 5 > 1, a low-temperature 
antiferromagnetic phase is predicted, T/v approaching the finite value (9), as 
J/W — > oo. This arises from the fact the 25-fold degenerate ground state set 
of the Kondo ion reacts with a Curie-like polarization to the transverse field 
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t x . The argument can be made a bit more formal by deriving the effective 
Hamiltonian which expresses the action of Hkn within the restricted Hilbert 
space spanned by the \S + t\ = S — 1/2 Kondo ground states. Expressed in 
the basis \S z ,t z ), the relevant single-site states are 

I i i 

5 -m-i/2 = v / 2ig+ X W^S -m\S - m, --) - y/m+\\S - m - 1, -)] (10) 

where m — 0, 1, 25"— 1. It is easy to check that in the Hilbert space spanned 
by these local basis states, Hkn acts like 

where the components of the usual spin operator r\ for \rj\ = S — 1/2 appear. 
Thus in the limit J/W — > oo, the Kondo necklace model becomes equivalent 
to an antiferromagnetic x — y model of the residual S — 1/2-spins, with an 
effective coupling W/ (2S+ 1) 2 . Note that this result is rigorous, being perfectly 
independent of the previous mean-field argument. 

It might seem unphysical that T N approaches a finite value as J/W — > oo, 
since it looks obvious that it should rather tend to zero. There is, however, 
no real contradiction. We have to remember that W is an intersite spin- 
flip matrix element of the conduction electrons, and is thus, in terms of the 
underlying Kondo lattice model (1), itself dependent on J. In particular, at 
large J, the r-spin-flip process involves breaking up the local Kondo ground 
states, and thus W oc B 2 / J is expected, where B is the conduction electron 
bandwidth. Hence (9) actually predicts T N — > as J/W — > oo. 

We now return to the mean-field solution, and discuss the low-temperature 
properties. Let us recall that Doniach's mean-field solution was actually done 
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for T = 0, being based on a site-factorized Ansatz for the ground state wave 
function. Our finite-temperature mean-field approximation should become 
equivalent to this at T = 0, since it relies on the same kind of decoupling. 
This becomes quite clear from a detailed study of the S — 1 case 0] where 
both formulations were implemented. However, for general S, we find it more 
convenient to use our present formulation. 

The solution of the eigenvalue problem of (6) reveals that the ground state 
is lying in the m = 2S subspace. The T = self-consistency equation has no 
closed-form solution. However, we can carry out series expansions, either for 
small, or large J/W. The large- J expansion up to fourth order gives 

1 2S- 1 16 S W 1QS (3 - 28S + 12S 2 ) W 2 

{T) ~22STl' { + (1 + 25) 3 T (1 + 2Sf I 2 ' 

1285 (1 - 285 + 136S 2 - 112S 3 + 16S 4 ) W 3 __ 
(1+2S) 9 ~J*~ 
64S (5 - 2905 + 3548S 2 - 129765 3 + 14192^ 4 - 4640S 5 + 320S 6 ) W A 

(1 + 2S) 12 J 4 " 1 (12) 

This is, roughly speaking, a power series in J/WS 2 , i.e., in the dimensionless 
effective coupling we can identify from (9). The overall factor 25 — 1 reflects 
that, as known from Doniach's early work 0, the magnetic order is suppressed 
by the formation of Kondo singlets for large J's if S = 1/2. On the other hand, 
no sign of any sudden change in the nature of the ground state is detected for 
S> 1. 

t x is perfectly valid as an order parameter but we rather prefer using S x , 
the transverse sublattice magnetization of the S-spins. Evaluating this needs 
no further calculation since, by construction 

(r x + S x ) = ±(S - 1/2) (13) 
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whereby it has to be observed that sg(S x ) = — sg(r x ). With our convention 
(3), the S-spin is "up" on the B sublattice. 
The fourth-order small-J expansion yields 

(r*) « - - - • + S(2S - 1} • ^ - -S(3 - 10S + 12S 2 ) ■ — A (14) 
N ' 2 2 W 2 2 W 3 8 V W 4 

S x for arbitrary J/W can be obtained numerically: examples are shown 
in Fig. 2. Generally, S — 1/2 < S x < S. The downward deviation from S 
is a measure of the strength of the (necessarily partial) Kondo compensation 
which depends not only on J/W but also (and rather drastically) on S. For 
large S, the S-spin density wave amplitude approaches the quasi-classical 
S — 1/2S, which is accompanied by an oppositely polarized T-spin wave of 
amplitude 1/2 — 1/2S. This picture of a frozen state yielded by the mean- 
field treatment apparently fails to do justice to the more dynamical internal 
structure of the strongly bound low-spin Kondo state, in terms of which the 
effective hamiltonian (11) is formulated. - To see where essential features 
may be missed let us note that in the mean-field treatment, "turning the spin 
of a site down" is achieved by acting with a linear combination of S~, and t~. 
However, the relevant operators are the ^-operators describing the composite 
Kondo-bound objects and in general, 7]~ will contain also more complicated 
components like S~S~t + , etc. 

3 Discussion and Summary 

Motivated by Doniach's study of the original S = 1/2 necklace model, 
we introduced the underscreened Kondo necklace model Hkn given in (2) 
for studying the situation where conduction electrons in a single screening 
channel can only partially compensate a lattice of localized spins S > 1/2. 



10 



The necklace model incorporates only one aspect of the full Kondo lattice 
problem: the competition between Kondo compensation, and intersite spin 
exchange, and is thus meant to model magnetic ordering in the case when the 
(non-degenerate) conduction band is half-filled, and the system is expected 
to have an insulating ground state. Though the denomination "necklace" is 
suggestive of Id, the above-described situation is present in all dimensions, 
and we have in mind really three-dimensional systems. Within a mean-field 
approximation, we found that the low-temperature phase is a two-sublattice 
x — y antiferromagnet for arbitrary strength of the Kondo coupling J. It also 
follows that the excitation spectrum is gapless, in the same sense as the ordered 
state of the original necklace model was found to be. 

Obviously, the Kondo necklace models do not have real "Kondo physics" in 
them, in the sense that low-energy electron-hole excitations, and with them 
the possibility of the characteristic non-analyticities at weak coupling, are left 
out. Kondo divergencies are expected to be most important when the (in the 
lattice case, collective) Kondo effect hinders, or at least greatly suppresses, 
ordering. Such is not the case with the underscreened necklace models where 
we find well-developed antiferromagnet ism. It can be argued |15| that in the 
cases of clearcut ordering, the eventual inclusion of true Kondo physics would 
shift the phase boundaries, but would not modify the overall appearance of 
the mean-field phase diagram. Our results should be considered in this spirit. 

In addition to missing out on the charge degrees of freedom, the necklace 
models differ in another significant way from the Kondo lattice models: they 
turn out to describe x — y type magnets rather than isotropic ones. This was 
already implicit in the fact that Doniach's || mean field trial state turned out 
to be polarized along the spin-x axis, but some subsequent studies |12| seemed 
to hint that the Kondo necklace may be, after all, in the same universality 
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class as the isotropic Kondo chain. This question is bound up with that of 
the existence of a spin gap in the spectrum. There seems to be no doubt that 
the spin excitation spectrum of the Id Kondo lattice is gapful ||13||. The recent 



numerical finding |L0j of a phase transition between a gapless and a gapful state 
in the ground state of the S = 1/2 Kondo necklace underlines two points: a) 
that the anisotropic form of the r-term in (2) has essential consequences, and 
b) that the mean field treatment is a useful guide to finding out whether there 
is such a phase transition. 

In the underscreened necklace models, the x — y character becomes even 
more emphatic. We derived the exact form (11) of the effective hamiltonian 
governing the large- J behaviour, and found it to be a pure antiferromagnetic 
x — y coupling between the residual composite spins. - From this it also 
follows that - in contrast to the expected behaviour of the Id underscreened 
Kondo lattice models - the Kondo necklaces show no Haldane phenomenon, 
i.e., the nature of the ground state is not expected to alternate between integer 
and half-integer values of S. This should be clearest at J/W — > oo when (11) 
can be used. In the ground state phase diagram of a larger class of spin models 



I6fl , the pure antiferromagnetic x — y model is lying on the boundary between 



the Haldane, and x — y phases, and is supposed to do nothing exotic. 

The results for the ground state sublattice magnetization are given in Fig. 2. 
It is a question of great current interest |7|, |2j , how the /-spins get apparently 
divided into a screened, and an ordered part. The simplification brought by 
replacing the Kondo lattice model with the necklace model allows a simple 
treatment of this problem. As J/W is increased, the order parameter shows 
a perfectly smooth behaviour, gradually decreasing from its unscreened value 
at J = 0, towards the asymptotic value (9). 

To summarize, to gain some insight into the difficult problem posed by 
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the underscreened Kondo lattices J7J], we introduced the underscreened Kondo 
necklace models (2). These are in the same relationship to the Kondo lattices 
as Doniach's |§ original necklace model to the S = 1/2 Kondo lattice. The 
simplification has its price: charge degrees of freedom are not considered, and 
the spin-rotational invariance of the original Kondo lattice model is destroyed. 
By deriving the exact form (11) of the effective hamiltonian governing the 
behaviour of the lattice model at very large Kondo coupling, we found that 
it describes an x — y antiferromagnet of the residual spins. In an effective 
field treatment, the ground state is always antiferromagnetic, and the spec- 
trum gapless. The mean field phase diagram is shown in Fig. 1. Though 
obviously oversimplified, the model introduced by us has the merit of allowing 
the description of the competition between Kondo compensation and magnetic 
ordering for underscreened localized moments. 
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Figure Captions 



Fig. 1 Neel temperature T N in units of W versus J/ (J + W), with d = 1, 
for several values of S. Note that for the underscreened models (S > 1/2), 
T N /W tends to a finite value as J/W — > oo. 

Fig. 2 The T = value of the order parameter S x changes in the range 
[S — 1/2, S] as a function of J/W (d — 1 was taken). For the underscreened 
S > 1/2 models, no ground state phase transition is found. 
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